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Abstract  This  paper  presents  an  integrated  missile  guidance  and  control  law  based  on  adaptive 
fuzzy  sliding  mode  control.  The  integrated  model  is  formulated  as  a  block-strict-feedback  nonlinear 
system,  in  which  modeling  errors,  unmodeled  nonlinearities,  target  maneuvers,  etc.  are  viewed  as 
unknown  uncertainties.  The  adaptive  nonlinear  control  law  is  designed  based  on  backstepping 
and  sliding  mode  control  techniques.  An  adaptive  fuzzy  system  is  adopted  to  approximate  the 
coupling  nonlinear  functions  of  the  system,  and  for  the  uncertainties,  we  utilize  an  online-adaptive 
control  law  to  estimate  the  unknown  parameters.  The  stability  analysis  of  the  closed-loop  system  is 
also  conducted.  Simulation  results  show  that,  with  the  application  of  the  adaptive  fuzzy  sliding 
mode  control,  small  miss  distances  and  smooth  missile  trajectories  are  achieved,  and  the  system 
is  robust  against  system  uncertainties  and  external  disturbances. 

©  2014  Production  and  hosting  by  Elsevier  Ltd.  on  behalf  of  CSAA  &  BUAA. 


1.  Introduction 

Missile  guidance  and  control  systems  are  usually  designed  sep¬ 
arately  due  to  the  assumption  that  there  is  a  spectral  separa¬ 
tion  between  the  guidance  loop  and  the  control  loop.  Based 
on  this  paradigm,  a  number  of  past  missile  systems  which 
guarantee  outstanding  performance  have  been  designed.  How¬ 
ever,  it  can  be  argued  that  this  design  paradigm  cannot  fully 
exploit  synergistic  relationships  between  the  two  subsystems 
or  strictly  maintain  the  stability  of  the  overall  system.1  On 
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the  other  hand,  the  spectral  separation  assumption  may  be 
invalid,  especially  at  the  end-game  phase  of  the  interception.2 
Integrated  guidance  and  control  (IGC)  design  was  first  put  for¬ 
ward  in  Ref.3,  and  has  received  much  attention  in  recent 
years.4-8  It  was  shown  that  IGC  designs  have  the  potential 
to  enhance  missile  performance  by  viewing  the  two  subsystems 
as  an  integrated  system  and  accounting  for  the  coupling 
between  guidance  and  control  dynamics. 

Various  control  methods  have  been  adopted  in  IGC 
designs.  A  small-gain  theorem  based  IGC  law  was  designed 
in  Ref.1  for  missiles  steered  by  both  canard  and  tail  controls, 
and  the  stability  of  the  overall  system  could  be  guaranteed 
without  the  assumption  that  the  angle  between  line-of-sight 
(LOS)  and  missile  velocity  was  almost  invariable.  An  IGC 
law  based  on  adaptive  output  feedback  and  backstepping  tech¬ 
niques  was  designed  in  Ref.7  for  formation  flight,  which  was 
translated  into  better  transient  and  steady-state  range  tracking 
performance.  An  IGC  law  based  on  the  state-dependent 
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Riccati  equation  approach  for  a  moving-mass-actuated  missile 
was  designed  in  Ref.8,  and  miss  distances  which  were  much  less 
than  the  diameter  of  the  missile  were  achieved.  The  nonlinear 
optimal  control  technique,  the  Q-D  method,  was  employed  in 
Ref.9  to  design  an  IGC  law,  and  the  controller  did  not  require 
online  computation  of  the  state-dependent  Riccati  equation. 

Sliding  mode  control  (SMC)  is  another  typical  method  in 
IGC  designs.  SMC  is  known  to  be  an  efficient  control  tech¬ 
nique  applicable  to  a  wide  class  of  nonlinear  systems,  due  to 
its  insensitivity  to  model  uncertainties  and  external  distur¬ 
bances  after  reaching  the  sliding  phase.  SMC  has  been 
addressed  in  some  previous  studies  for  IGC  designs. 210-14 
Koren  et  al.2  chose  the  zero-effort  miss  distance  as  the  sliding 
variable.  A  robust  SMC  controller  was  then  designed  to  deal 
with  both  system  uncertainties  and  the  difference  between  non¬ 
linear  and  linear  design  methods.  Shima  et  al.10  defined  the 
same  sliding  surface  as  that  in  Ref.-.  Based  on  their  approach, 
small  distances  could  be  achieved  even  in  stringent  interception 
scenarios.  Hou  and  Duan1 1  proposed  an  IGC  scheme  for  hom¬ 
ing  missiles  against  ground  fixed  targets,  and  an  SMC-based 
adaptive  nonlinear  control  law  was  designed  to  guarantee  a 
missile  hit  a  target  accurately  with  a  desired  impact  attitude 
angle.  Based  on  the  assumption  that  each  of  the  three  channels 
of  an  IGC  model  can  be  independently  designed,  Yamasaki 
et  al.12  introduced  an  IGC  design  approach  for  a  path¬ 
following  uninhabited  aerial  vehicle.  Dong  et  al.13  developed 
a  robust  higher-order  sliding  mode  (HOSM)  based  IGC  law, 
in  which  the  IGC  design  problem  was  considered  to  be  equal 
to  the  stabilization  of  a  third  integral  chain  system.  Zhao 
et  al.14  proposed  a  SMC-based  nonlinear  IGC  strategy  which 
took  the  higher-order  dynamics  of  the  system  into  account. 

Although  SMC  has  been  widely  applied  to  IGC  designs, 
some  problems  still  exist.  Nearly  all  existing  approaches  are 
based  on  the  assumption  that  the  nonlinear  functions  in  an 
IGC  model  could  be  accurately  obtained.  In  practice,  such 
an  assumption  may  not  be  always  guaranteed.  In  this  paper, 
an  IGC  law  based  on  adaptive  fuzzy  sliding  mode  control  is 
firstly  presented.  The  developed  approach,  when  compared 
with  the  existing  results,  is  novel  in  that  the  IGC  law  can 
guarantee  high  performance  without  the  assumption  that  the 
coupling  nonlinear  functions  in  the  integrated  model  can  be 
accurately  obtained. 

2.  Model  derivation 

2.1.  Engagement  kinematics 

The  planar  engagement  geometry  is  depicted  in  Fig.  1,  where 
OXY  is  a  Cartesian  inertial  reference  frame,  and  M  and  T 
represent  the  missile  and  the  target,  respectively.  The  corre¬ 
sponding  equations  of  motion  between  the  missile  and  the 
target  are  as  follows:1 

R  =  Ft  cos {q  -  0T)  -  VM  cos {q  -  dM)  (la) 

Rq  =  —V x  sin(<7  —  0T)  +  VM  sin(^  —  0M)  (lb) 

where  R  is  the  relative  range,  q  is  the  LOS  angle,  0M  and  0T  are 
the  missile  and  target  flight  path  angles,  respectively,  and  VM 
and  Vt  are  the  missile  and  target  velocities,  respectively. 
Differentiating  Eq.  (lb)  followed  by  the  substitution  of 
Eq.  (la),  we  get 


Fig.  1  Planar  engagement  geometry. 


Rq  +  2 Rq  =  -VT  sin (q  -  0T)  +  VM  sin(?  -  0M)  +  VT0T 

x  cos (q  -  0T)  -  Vm0m  cos (q  -  0M)  (2) 

Assume  that  VM  =  VT  =  o,  and  define  Vq  =  Rq,  aT  =  Vy  Oj, 
and  =  Pm  Om  •  Eq.  (2)  can  be  rewritten  as 

Vq  =  -  ^  Vq  +  aT  cos (q  -  0T)  -  aM  cos (q  -  0M)  (3) 

where  aM  and  aT  are  the  missile  and  target  accelerations, 
respectively. 

2.2.  Missile  dynamics 

The  planar  missile  dynamics  are  given  by15 


«  =  m  T’m  sin  *  —  L  +  mg  cos  0M)  +  wz  (4) 

Jzwz  =  M0  +  MSz, 5_.  (5) 

*  =  (6) 

«  =  (7) 


where  a  is  the  angle  of  attack,  m  is  the  missile  mass,  TM  is 
the  thrust  of  the  missile,  L  is  the  lift  force,  o>z  is  the  pitch 
rate,  Jz  is  the  moment  of  inertia  about  z-axis,  Sz  is  the 
deflection  angle  for  pitch  control,  0  is  the  pitch  angle,  Ms,  is 
the  control  contribution  to  the  angular  acceleration,  and 
M0  =  M0(ot,  Ma,  h,  VM,  coz)  represents  the  angular  accelera¬ 
tion  contributions  from  all  other  sources  such  as  the  angle  of 
attack  a,  the  Mach  number  Ma,  the  height  h,  and  so  on.  M0 
is  often  approximated  as  follows:16 

M0  =  Ma<x.  +  Mm,mz  (8) 

where  Ma  and  Mm,  are  the  angular  acceleration  contributions 
from  the  angle  of  attack  and  pitch  rate,  respectively. 

The  lift  force  (L)  and  relative  parameters  (Mz,Mm,,Ms,) 
are  as  follows: 

IL  =  51.?>Qs(cyU.  +  c^<5z) 

Ma  =  51 3Qslm*rj. 

M  _  Qstm f  W 

Ms,  =  51.3Qslmz‘ 

where  Q  is  the  dynamic  pressure,  s  is  the  aerodynamic  refer¬ 
ence  area,  /  is  the  reference  length,  c“  and  c*  are  the  lift  force 
derivatives  with  respect  to  a.  and  <5Z,  respectively,  and  mJ2,  m'"z, 
and  mf  are  the  pitch  moment  derivatives  with  respect  to  a,  a>z, 
and  bz,  respectively. 
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2.3.  Integrated  model 


According  to  the  above  analysis  and  with  the  assumption  that 
during  the  end-game  the  missile  has  no  thrust  and  the  drag 
force  is  negligible,  the  integrated  model  can  be  written  as 


V*=-RV*  +  - 
57.3  Qsc“ 


- 51.3Qsc ; 


g cos  9m 

mVM  ”  '  VM  + 
_  57 3Qslm“  _  ,  57 3Qsfm 
Jz 

573Qslm\ } 


JzV M 


Jz 


—  Sz  +  d(il2 


where  dVq ,  dy,  and  da>2  are  the  unknown  bounded  uncertainties 
(modeling  errors  caused  by  the  assumptions,  unmodeled  non- 
linearities,  target  maneuvers,  etc.). 

V„ 

Define  xi  =  ,  x2  =  a,  x3  =  mz,  u  =  dz,  and 

„„„  ,  *  -57 3Qsc,Jm 
.  57 3Qslm°*  .  y 

b  = - — —  .  The  integrated  model  can  be  rewritten  into 

fit  =fi(xi)  +  x2  +  di 

X2  =  f2(x2)  +  X3  +  d2  (11) 

x3=f3(x2,x3)  +  bu  +  d3 
where 


J\(xi)  =  ~kxi  + 
-57.3  Qsc; 


R  gcosdMcos(q  —  9m) 


— 57 3Qsty/m 


f3(x2,x3)  = 
di  = 


mV m  z  '  VM 

57 3Qsbn*  Qsfrrt 

Jz  ** +  JzW y 
d2  =  d. 


57 3Qsc“/m  ’ 
d3  =  dm„\di\  ^dmax  (i  =  1,2,3) 


estimated  value  f\(x\).  Due  to  the  universal  approximation 
ability  of  an  adaptive  fuzzy  system,  it  will  be  used  here  to 
approximate  the  uncertain  nonlinear  function  f\{x\).  For  a 
continuous  function  f(x)  defined  in  a  close  set  U  and  any 
precision  g,  there  must  exist  a  fuzzy  logic  system  F(x)  that 
satisfies18,19 

sup|F(jc)-/(*)|<e  (14) 

xeu 

where  F(x)  is  composed  of  product  inference,  singleton  fuzzifi¬ 
cation,  center-average  defuzzification,  and  Gauss  membership 
functions. 


(1)  Product  inference 

\t0)  =  /y^O/y^)  •  ■  ■  /y„(*»)  (15) 

where  x  =  [xi,x2, . . .  ,x„]T  and  f  are  the  fuzzy  logic  system’s 
input  and  output,  respectively.  Fuzzy  sets  /  and 
ft{i  =1,2 ,...,«)  are  associated  with  the  fuzzy  membership 
functions  /(}(y/)  and  /y  (xf)  (i  =  1,2 ,...,«),  respectively. 

(2)  Singleton  fuzzification 

The  rules  to  map  a  crisp  point  x  =  [xi,x2, . . .  ,x„]T  G  U into  a 
fuzzy  set  Ax  in  U  are:  Ax  is  a  fuzzy  singleton  with  support  x, 
i.e.,  n a  (V)  =  1  for  x!  =  x  and  Ha(x)  =  0  for  all  other  x'  €  U 
with  x1  &  x. 

(3)  Center-average  defuzzification 


EUf-m 

y  TzUm 


(16) 


(4)  Gauss  membership  function 


/y(*i)  =  exp  \y  j  (17) 

where  x\  is  the  center  of  the  Gauss  curve  and  <rf  is  the  width  of 
the  Gauss  curve. 


By  defining  the  fuzzy  base  function 


With  simplicity  consideration,  the  IGC  model  (11)  is  in  a 
strict-feedback  form,17  and  the  backstepping  technique  is  suit¬ 
able  for  dealing  with  such  cascade  systems.  Moreover,  SMC  is 
known  to  be  an  efficient  control  technique  to  overcome  model 
uncertainties.  The  main  observation,  in  this  paper,  is  to 
combine  these  two  techniques  to  design  a  high-performance 
IGC  law. 

3.  IGC  law  design 

The  task  of  the  IGC  law  design  is  to  find  a  suitable  control  u  to 
make  sure  the  missile  hit  the  target,  and  during  the  process  the 
missile  attitude  is  stable.  The  IGC  law  design  procedure  in  this 
paper  is  divided  into  three  steps. 

Step  1.  Based  on  the  intuition  that  zeroing  q  will  lead  to 
interception,  we  choose  the  first  sliding  surface  as 

Si  =  xi  (12) 

The  time  derivative  of  ,S|  is 

si  =*i  =fi(xi)  +  x2  +  di  (13) 

In  the  control  system,  we  cannot  measure  the  exact  value  of 
the  nonlinear  function  /i(xi),  so  it  will  be  replaced  by  its 


the  fuzzy  system  can  be  written  into  the  following  equivalent 


W  f  „  s  =  ep(x) 


(19) 


where  0  =  \y\f, . . .  ,r]T,P{x)  =  [Pi(x),  P2{x), ...,  Pm(x)]T. 

Details  about  this  kind  of  adaptive  fuzzy  systems  can  be 
found  in  Refs.18,19.  In  this  paper,  we  only  need  to  construct 
the  approximate  Gauss  membership  functions  of  the  system 
states.  In  practice,  the  membership  functions  will  be  deter¬ 
mined  by  an  iterative  procedure  according  to  the  computa¬ 
tional  results. 

Assuming  an  adaptive  fuzzy  system 


Ax)  =  0jP/x) 


(20) 


where  9/  is  an  unknown  weight  vector  that  needs  online  regu¬ 
lation  and  Pj{x)  is  a  Gauss  function  as  well  as  a  fuzzy  base 
function. 
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During  the  approximation  process,  we  can  get  the  optimum 
value 


fix)  =  0}JPj{x) 

(21) 

where  0*f  is  the  optimum  weight  vector  defined  as 

0}  =  argmin  [sup|/(x)  -/Ml] 

(22) 

Yxtu  \ 

where  /  is  the  estimate  of  /.  Note  that  we  will  us 

se  the  corre- 

sponding  hat  “A”  to  denote  the  estimate  value  henceforth. 

As  shown  in  Eq.  (14),  there  exists  an  approximation  error 
that  satisfies 

/(*)=/(*) +  e/  (23) 

where  |e/|  <  rff,  in  which  r\*f  is  the  upper  boundary  of  the 
approximation  error.  However,  the  optimum  parameter  0)  is 
also  substituted  with  the  estimated  value 

m  =  OjPf{x)  (24) 

From  the  above  analysis,  we  have 


where 

!02  =  rj\s2P2{x2) 

h  =  4\M  (32) 

^2, max  =  rfi\s2\ 
with  17?  >0  (;'  =  1,2,3). 

In  Eq.  (31),  x2c  can  be  obtained  from  a  first-order  low-band 
filter.  With  the  filter,  the  “computation  explosion”  problem  in 
traditional  backstepping  designs  can  be  avoided.20 

Step  3.  The  third  sliding  surface  is  given  as 
ft  =  -ft  -  x3c  (33) 

The  time  derivative  of  ft  is 

ft  =  x3  -  X3c  =  /3(x2,  x3)  +  bu  +  d3  —  x3c  (34) 

In  this  step,  we  obtain  the  actual  control  signal 
u  =  —b~l(s2  +  k3s3  +  0jP3(x2,  x3)  +  i3sgn(ft) 

+  4,maxSgn(ft)  -  X3c)  (35) 

where 


/(*)  =f(X)  +  Bf=Kx)  +/j»  +8/=  OJP/X)  +  OjP/x)  +  8f 

(25) 

where  f(x)  =  f{x)  — f(x )  and  Of  =  0"f  —  Of  are  the  approxima¬ 
tion  error  and  the  weight-value  approximation  error, 
respectively. 

Therefore,  with  the  result  in  Eq.  (25),  Eq.  (13)  can  be 
rewritten  as 

ft  =fi(xi)  +/i(xi)  +  El  +  x2  +  d\ 

=  OjPi  (. *1 )  +  OjPi  (xi)  +Bi+x2  +  dl  (26) 

We  treat  x2  as  the  virtual  command,  which  can  be  designed 


X2c  =  — kift  -  b]P\  (xi)  -  ftSgn(ft)  -  di,maxsgn(ft)  (27) 

For  the  uncertain  system  parameters  0i,ei,  and  d\ jmax  are 
unknown,  the  following  online-adaptive  control  law  is 
proposed: 

f  k  =  flUlPl(Xi) 

(28) 

l^l,max  =  »/3|ft| 

where  r/j  >  0  (i  =  1, 2, 3). 

Step  2.  The  second  sliding  surface  is  given  as 
ft  =  x2-x2c  (29) 

The  time  derivative  of  ft  is 

s2  =  x2-x2c=f2{x2)  +  x3  +  d2-x2c  (30) 

In  this  step,  we  choose  x3  as  the  virtual  command,  and  the 
derivation  steps  are  quite  similar  to  those  in  Step  1.  For  the 
sake  of  brevity,  we  give  x3c  directly: 

-v3c  —  -  ft  -  k2s2  -  blP2(x2,x3)  -  s2sgn(ft) 

-  4maxSgn(ft)  +  ftc  (31) 


{b3  =  t]]s3p3(x2,x3) 

b-,lh\  (36) 

^3,max  =  ^|ft| 

b~l  —  r\\s3bu 

with  >/?  >  0  (i  =  1,2, 3, 4). 

Finally,  we  state  the  complete  control  law  as  follows: 

’  ft  =  Xi 

X2C  =  —klSi  -  0]Pi(Xi)  -  ftSgn(ft)  -  (J^maxSgn^!) 
ft  =  x2  -  x2c 

*3c  =  --Si  -  k2s2  -  0jP2(x2,  x3)  -  e2sgn(ft) 

-^2,maxSgn(ft)  +  X2c 

ft  =  X3  -  X3c 

u  =  —b~l(s2  +  k3ft  +  0jP3{x2,  x3)  +  83sgn(ft) 

+<?3,maxSgn(i3)  -  X3c) 

The  online-adaptive  laws  of  the  unknown  parameters  are 
defined  in  Eqs.  (28),  (32),  and  (36). 

4.  Stability  analysis 

Define  0,  =  0,  —  0,-, s,-  =  e;  —  ef,  dim.M  =  rfimax  —  c4max  [i  =  1 , 2, 3), 
and  b~1  =  b~l—b~l.  Substitute  Eqs.  (27)  and  (29)  into 
Eq.  (13),  and  get 

ft  =fi(xi)  -  Aqft  -  0[Pi(x,)  -  Eisgn(ft)  -  ^i,maxsgn(ii)  +  d\ 

+  ft  =  0[Pi(xi)  +  si  -  kift  -  Eisgn(^i) 

-  <7i,maxSgn(ft)  +  d\  +  ft  (38) 

Substitute  Eqs.  (31)  and  (33)  into  Eq.  (30),  and  get 

ft  =f%{x2)  -  ft  -  k2S2  -  0[P2{x2)  -  £2Sgn(ft)  -  rf2,maxSgn(ft) 

+  X2c  +  d2-  x2o  +  ft  =  02  P2(x2)  +  e2  -  ft  -  k2s2 
-  E2sgn(ft)  -  <72,maxsgn(ft)  +  d2  +  ft 


(39) 
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Substitute  Eq.  (35)  into  Eq.  (34),  and  get 
Is  =/3  ( x2  ,  x3 )  -  bb~ 1  (s2  +  k3s3  +  03  P3{x2,  x3 )  +  e3  sgn(j3 ) 

+  d3,maxSgn(j3)  -  X3c)  +  rf3  -  X3c  =/3(x 2,X3) 

-  {s2  +  k3s3  +  03  P3(x2,x3)  +  83sgn(s3)  +  d3,maxsgn(i3)  -  x3c) 

+  bb~'  (s2  +  k3 S3  +  0]P3(x2,X3)  +  83Sgn(s3)  +  d3,maxSgn(.S3)  -  x3c) 
+  d3-  x3c  =  83  +  0]P3(x2lx3) -s2-  k3s3  -  e3sgn(s3) 

-  d3,maxSgn(i3)  +bb~\s2  +k3s3  +  0jP3{x2,x3)  +  83sgn(.s3) 

+  d3,maxSgn(i3)  - x3c)  +  d3  =  e3  +  03 P3(x2,x3)  -s2 

-  k3s3  -  e3sgn(j3)  -  d3,maxsgn(s3)  -  bb~'bu+d3  (40) 


Theorem  1.  Considering  the  nonlinear  system  Eq.  ( 1 1)  with 
bounded  uncertainties,  if  adaptive  fuzzy  systems  are  utilized  to 
approximate  the  uncertain  functions  /i(xi),/2(x2),  and  f3(x 2,  x3), 
and  in  the  SMC  control  law  Eq.  (37),  the  adaptive  parameters 
are  adjusted  by  the  system  online  adaptive  control  laws  shown  in 
Eqs.  (28),  (32),  and  (36),  then  the  system  sliding  modes  are 
asymptotically  accessible  and  the  closed-loop  system  is  asymp¬ 
totically  stable. 

Before  the  proof,  we  need  the  following  lemma. 

Lemma  1.  [Babalat’s  Lemma]  If  fit)  is  a  uniformly  continuous 
function  and  lim^oo  f0f[ T)dt  exists,  then  f[t)  converges  to  zero 
asymptotically. 


Proof  of  Theorem  1.  Choose  the  quasi  Lyapunov  function  a 

lA,  1  1^1, 

>'-2&  +  2^,'e‘+2^vf 


Differentiate  V  along  the  track  of  the  system  (11),  and  obtain 


Substitute  Eqs.  (38)-(40)  into  Eq.  (42),  and  get 

V=Si{0\P\{x\)  +  ei  -A:isi-gisgn(ff1)-rfiiniaxsgn(ji)  +  rfi  +s2) 
+  i2  (oJP2{x2,x3)  +  e2  -  ji  -  k2s2  -  £2sgn(52) 

-4,maxSgn(52)  +  +  S3) 

+53  (o]P3{x2,x3)+e3  -s2-k3s3-  E3sgn(53) 

-  d3,max  sgn(s3 )  +  bb~ 1  bu  + 

—  E  —  0]  0;  —  E  — -2,-ef  —  E  ~ ~  f  max  di  max  +  —  6“ '  6“ '  < 
4^0!  tf02  ifhS  "  '  04 

-  kiS  \  +  Ji  0\PX  (x3  )  +  El  I  Ji  I  -  6,  |  Ji  I  +  di  ,max  1 5]  I  -  di  ,max  |si  I 

—  k2g  +  S202  P 2(x2,X3)  +  £2|s2|  —  82|j2|  +  rf2jmax|52|  —  ^.maxl^l 
- k3s]  +  s30]P3(x2,x3)  +  S3 |i3 1  -  §3 153|  +  (/3,max|53 1  -  d3tnax\s3 1 

—  S3bb~lbu  —  E  —  0] hi  -  E -fe&i  —  E -T^imax^imax 

+  bnb-'b-' = -jyj + 


+s202P2(x2,x3)  +  530jP3(x2,x3)  —s3bb  lbu-J2-Oj0i 
-  E  1  8,8, '  -  E  '  f, max  f, max  +  ^  b~' b~l  =  -'S^kifi 

ja'fi  ffri  ’  ’  04  tf 

+  E<W  (fct -^Eax^)  -  bb-1  s3bu -  (43) 

Associating  with  Eqs.  (28),  (32),  and  (36),  we  have 

r<  _£***<  °  (44) 


Thus  5,-,  0i.  s/,  rf,;raax  (i  =  1,2,3),  and  b  1  are  all  bounded. 

Define  k  =  vnin{k\,k2,k3}  and  s  =  [.s’ i  ,.y2,  s3]T.  Eq.  (44)  can 
be  rewritten  as 


(45) 

Integration  of  Eq.  (40)  from  t  =  0  to  t  — »  oo  reveals  that 


jT°  *||*(t) ||*dr  <  -  J™  V(x)dT  =  E(0)  -  K(oo)  <  +oo  (46) 

Applying  Lemma  1  to  Eq.  (43),  we  get  that,  while  t  —>  0, 
||*||2  -*■  0,  i.e„  St  0  (i  m  1,2,3). 

Therefore,  it  can  be  concluded  that  the  system  sliding 
modes  are  asymptotically  accessible  and  the  closed-loop 
system  is  asymptotically  stable,  and  this  completes  the 
proof.  □ 


5.  Simulation  results 


This  section  presents  simulation  results  of  the  proposed  IGC 
law  on  a  numerical  example  introduced  in  Ref.:l.  In  this 
simulation  study,  the  constant  missile  speed  is  assumed  to  be 
VM  =  3.5 Ma.  The  initial  missile  attitude  and  control  fins 
are  a  (0)  =  0°  and  <5,(0)  =  0°.  The  initial  missile  flight  path 
angle  is  0M(0)  =  0°.  The  constant  target  speed  is  assumed  to 
be  Ft  =  900  m/s.  The  initial  target  flight  path  angle  is 
0-i  (0)  =  10°.  The  initial  missile  position  coordinate  is  (0,16)  km. 
The  initial  target  position  coordinate  is  (1,16.4)  km.  The  mis¬ 
sile  model  parameters  are  as  follows: 


The  actuator  dynamic  is  approximated  as  a  first-order  time 
delay  system  with  a  time  constant  of  0.01s,  and  the  control  sur¬ 
face  deflection  limit  is  |<5Z|  <  30°. 

The  controller  parameters  are  given  as  ki  =  0.8,  k2  =  1.0, 
k3  =  1.0,  £-'(0)  =  -0.1,  0,(0)  =  [0,0, 0]T,  6,(0)  =  0,  and 
fmax(0)  =0  (i—  1,2,3).  Assuming  that  |xi|  <  2 it,  \x2\  < 
and  |x3|  ^  7t,  we  choose  the  fuzzy  base  functions  as  follows: 


Fig.  2  Simulation  results  of  the  proposed  IGC  law  in  Case  1 . 
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Fig.  3  Simulation  results  of  the  proposed  IGC  law  in  Case  2. 
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ftJ.«p(-(a)!),fcl«„P(-(s+25) ) 
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We  evaluate  the  proposed  IGC  law  in  the  following  two 
cases: 

Case  1:  Suppose  d,  =  0  (i  =  1,2,3),  which  means  the  target 
does  not  maneuver  and  there  are  no  external 
disturbances. 

Case  2:  Assume  that  the  target  escapes  with  an  acceleration 
of  c?T  =  3gsin(7it/3),  the  missile  aerodynamics  coeffi¬ 
cients  vary  +25%,  and  external  disturbances 
d.2  =  0.1  sin  t  and  d3  =  0.2sin  t. 

Simulation  results  of  Case  1  and  Case  2  are  depicted  in 
Figs.  2  and  3,  respectively.  From  Figs.  2(a)-(c)  and 
Figs.  3(a)-(c),  we  can  see  that  the  missile  flight  trajectories 
are  smooth,  the  miss  distances  we  get  are  quite  small 
(0.0749  m  and  0.1334  m  in  Case  1  and  Case  2,  respectively), 
and  Vq  rapidly  converges  to  a  small  neighborhood  of  zero. 
Figs.  2(d)-(e)  and  Figs.  3(d)-(e)  depict  the  responses  of  missile 
dynamics.  Figs.  2(f)  and  3(f)  show  the  histories  of  the  control 
deflection.  It  is  observed  that  the  maximum  control  effort  is  less 
than  20°.  Finally,  histories  of  the  sliding  surfaces  are  shown  in 
Figs.  2(g)-(h)  and  Figs.  3(g)-(h).  The  sliding  surfaces  converge 
to  nearly  zero  rapidly,  which  guarantees  the  interception. 

6.  Conclusions 

(1)  With  some  rational  assumptions,  the  IGC  model  can  be 
built  in  a  strict-feedback  form,  and  by  adopting  the 
backstepping  technique,  the  stability  of  the  entire  system 
states  is  guaranteed. 

(2)  The  adaptive  fuzzy  system  is  effective  to  approximate 
the  coupling  nonlinear  functions  and  the  online-adaptive 
control  law  is  suitable  to  estimate  the  unknown  param¬ 
eters  in  the  integrated  system. 

(3)  Simulation  results  confirm  the  effectiveness  of  the  proposed 
method  on  dealing  with  missile  aerodynamics  coefficients 
varying,  target  maneuver,  and  external  disturbances. 
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